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We recently proposed a scheme for the creation of coherent superpositions of vortex states in Bose-Einstein 
condensates (BEC) using orbital angular momentum (OAM) states of light [Phys. Rev. Lett. 95, 173601 
(2005)]. Here we discuss further technical details of the proposal, provide alternative, time-reversal-symmetric 
scheme for transfer of a superposition of OAM states of light to the BEC via a procedure analogous to the tradi- 
tional STimulated Raman Adiabatic Passage (STIRAP) technique, and discuss an alternative trap configuration 
conducive for sustaining large charge vortices. Superpositions of OAM states of light, created using experi- 
mental techniques, can be transfered to an initially nonrotating BEC via a specially devised Raman coupling 
scheme. The techniques proposed here open up avenues to study coherent interaction of OAM states of light 
with matter. The study could also be employed for performing various quantum information processing tasks 
with OAM states of light — including a memory for a quantum state of the initial superposition. 
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I. INTRODUCTION 

The interaction of light, carrying quantized orbital angu- 
lar momentum (OAM) JjJ], with atomic or molecular matter 
is of considerable research interest i,iH. This quest is in- 
creasingly becoming feasible due to the progress that has been 
made in creation iOU^I 01, manipulation |8], detection ifl floll . 
and application [ 1 1 ] of the orbital angular momentum (OAM) 
states of light. Nevertheless, due to the size-mismatch be- 
tween the atoms and the spatial features of the phase struc- 
tures in the OAM beams of light, it is difficult to couple the 
OAM degrees of freedom of light to the internal states of 
single atoms. As a result, some of the activity in this area 
is concentrated in the excitation of vortices in Bose-Einstein 
Condensates lll2lll3lll4ll . Bose-Einstein Condensates (BEC) 
are macroscopic coherent objects that are obtained by trap- 
ping atoms in particular internal states and cooling them fur- 
ther so that all the atoms have the same motional and internal 
states. Excitation of vortices in BECs have been traditionally 
achieved via stirring of the BEC cloud with a laser beam 1151 . 
These vortex states are fairly stable and could be candidates 
for qubits in Quantum Information [ 16], if appropriate means 
to manipulate them are developed. 

There exist several proposals to transfer pure OAM states 
of light to the BECs 1171 [1911 : there is also a proposal to use 
BECs for storage of the OAM state of light J20[]. This transfer 
of OAM of light to BEC is possible due to the coherent nature 
of the BEC and because the size of the BEC cloud is about 
the same as the distance over which the intensity variations 
occur in the beam of light carrying OAM. We recently pro- 
posed a scheme for creation of superposition of two counter- 
rotating vortex states (with same absolute value of the charge) 
in BECs via coupling of the OAM degrees of freedom of inci- 
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dent beam to the center-of-mass motion of the atoms Ir2lll . To 
illustrate, transfer of a pure OAM state of light to a BEC cloud 
is an experimental reality [22]; nevertheless the experimental 
system involves different linear momentum for the initial and 
final states of the BEC as opposed to different internal states 
in the scheme discussed in Ref. I2H1 and in this article. Fur- 
thermore, the experiment involves transfer of only pure OAM 
states to the BEC as opposed to a OAM superposition as dis- 
cussed here. It is important to note that coherent transfer of a 
superposition state is the first check necessary to test the suit- 
ability of these macroscopic objects for quantum information. 
The results of our study reported here demonstrate the possi- 
bility of using vortex states in BEC as a qubit. 

The counter-rotating vortex superposition generated via our 
scheme is analogous to the counter-rotating persistent cur- 
rents in superconducting circuits |23[ |24[ |25|] . These counter- 
rotating persistent currents in the superconducting systems 
have already been recognized as candidates for qubits for 
quantum information processing. Thus the application of the 
study presented here to the area of quantum information is 
clear. Further applications of our scheme could be as a quan- 
tum memory of such superposition states. One can envision 
such a memory device to be useful in the quantum communi- 
cation networks using OAM states of light. (See Ref. B26I1 and 
references therein.) OAM states of light offer a higher dimen- 
sional Hilbert space to obtain extra security and dense cod- 
ing of quantum information. The memory application would 
deem it necessary to have a time-reversal symmetric proce- 
dure to write the superposition in OAM states of light to BEC 
so that the superposition can be extracted at a later stage. 
Here we discuss further technical details of our earlier pro- 
posal l2lll and also provide a STIRAP-like time-reversal sym- 
metric mechanism l27ll for transfer of special OAM states of 
light to BEC. Furthermore, we discuss alternative trapping po- 
tential, which can support toroid-shaped condensate cloud, in 
order to increase the stability of the generated vortex superpo- 
sition. 

The article is organized as follows. In Sec. [II] we discuss 
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generation of arbitrary superposition of two counter-rotating 
OAM components with the same absolute value of charge. 
This section also discusses a scheme for generation of arbi- 
trary superposition of two arbitrarily charged OAM states of 
light. Next, in Sec. [Til] we discuss two independent methods 
of transferring the OAM superposition from the light beam 
to the vortex states of BEC in complete detail. We also con- 
sider an alternative trap configuration, the so-called mexican- 
hat trap to offer an example of trap that can sustain large 
charge vortices. In fact, we provide a very general set of equa- 
tions to study the transfer of OAM superposition from light to 
BEC for any kind of trapping potential. In Sec.[IV]we discuss 
the detection of the BEC vortex superposition and present our 
conclusions in Sec.[V] Some of the mathematical details are 
presented in Appendices A and B. 



II. SUPERPOSITION OF OPTICAL VORTICES 

To recollect, the OAM states of light are different in their 
phase characteristics from the conventional gaussian light 
beams. They have azimuthal phase structure characterized by 
the quantized orbital angular momentum carried by each pho- 
ton in the beam. In the simplest case of one unit of orbital 
angular momentum, i.e. h, the light beam has a cork-screw 
type phase front. Thus, the phase continuously varies from 
to 2n along the azimuthal coordinate and there is a jump from 
2tz to along some radial line. The location of the radial line 
along which the phase jump (or phase discontuinity) occurs 
is time-dependent and it rotates continuously with time. The 
sense of rotation governs the sign of the the orbital angular 
momentum. Furthermore, an OAM state with angular mo- 
mentum th has \£\ number of azimuthal phase jumps across a 
cut taken in the beam path. It has to be noted that this is an 
instantaneous description of the beam phase across a cut taken 
in its path. The whole structure is time dependent and rotates 
along the beam axis in a clockwise or counter-clockwise di- 
rection as the beam propagates. The sign of £ corresponds 
to the sense of rotations of the phase fronts around the beam 
axis. In principle, there is no upper limit to the angular mo- 
mentum value that can be imparted to a light beam. Therefore, 
OAM states of light themselves are very good candidates for 
quantum information processing, mainly in the area of quan- 
tum cryptography, as they offer an essentially infinite Hilbert 
space to work with. This could result in increasing the secu- 
rity of the quantum-key distribution protocols tremendously. 

Monochromatic beams with azimuthal phase dependence 
of the type exp(i£0), of which Laguerre-Gaussian (LG) laser 
modes are an example, have well-defined angular momentum 
of th per photon. The normalized Laguerre Gaussian mode at 
the beam waist (z = 0) and beam size wo at the waist is given 
in cylindrical coordinates (p, <j>,z) by 

y n{\£\+p)\w y w J V w o ) 
exp(-p 2 /wo)exp(i€0), (1) 



where lk (p ) are the associated Laguerre polynomials, 

4V) =£(-D" , ^ „ , P"- (2) 
m~0 {p-m)\(\£\+m)\m\ 

Here, p is the number of non-axial radial nodes, and the index 
I, the winding number, describes the helical structure of the 
wave front and the number of times the phase jumps occur as 
one goes around the beam path in the azimuthal direction. In 
general, any mode function \j/(p,<j),z) can be expanded into 
LG modes as 

WM,z)= E X>^ LG p(P,<M (3) 

£=-<*> p=0 

For further discussion we consider only pure LG modes with 
a definite charge I and p = 0, we denote such a state of the 
light field by \l) such that (r\£) = LG$(p,0). Thus it can be 
easily seen that the states \+£) and \ (with £ being a whole 
number) differ only in the sense of the winding of the phase; 
either clockwise or counter-clockwise. 

Now we discuss the first step of the proposal: creation of 
an arbitrary superposition of two OAM states of light of the 
kind (a \£) + ]8 \£' ))/V« 2 + P 2 - First we discuss generation 
of a special kind of superposition state: (a \£) + j3 |— £)), with 
a 2 + j3 2 = 1. Here the two components of the state have op- 
posite sense of rotation. It is well known that creation of a 
very general multicomponent superposition of OAM states of 
the kind Y.£ c e is a fairly straightforward procedure by us- 
ing holographs or phase plates. Furthermore, a sorter of these 
OAM states has also been demonstrated 111 Oil , which distin- 
guishes between different orbital angular momentum compo- 
nents. Thus, by using a mixed OAM-state generator in con- 
junction with an OAM-sorter one can easily obtain a pure 
OAM state with arbitrary \£). 

Next, we describe a technique that could be used to trans- 
form a \+£) state in to a \—t) state with £ being a whole 
number. We illustrate this schematically in Fig. [TJ A com- 
mon representation of the OAM states of light is through the 
Laguerre-Gaussian beams which have circular cross-sections. 
We consider two rays (red and blue) as seen in Fig. [TJthat lie 
diametrically opposite to each other on the beam profile. As 
can be easily shown from the ray-diagram these rays would 
interchange their places in the beam profile l28l l29ll . Extrap- 
olating this to the whole beam it can be seen that passing a 
OAM state of light through a dove prism amounts to a anti- 
symmetrization of its azimuthal phase structure thus giving 
rise to a |— I) state. It is also important to note, as pointed out 
by Padgett et al. [29], that this rotation of the OAM state does 
not affect the polarization state of the light. To obtain com- 
plete conversion in the handedness of the input OAM, via the 
Dove prism, further care may be necessary as pointed out in 
Ref.il. 

Thus, starting with a \+£) state of light and passing it 
through a Mach-Zehnder type configuration shown in Fig. [2] 
we obtain ( | I) + \ — I) ) / y/2 at one of the output ports of the in- 
terferometer. By choosing the first beam splitter in the Mach- 
Zehnder interferometer as an a//3 beam splitter and the sec- 
ond one is a 50/50 beam splitter one can generate a general 
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FIG. 1 : Dove prism as a sign shifter for the OAM states of light. The 
dove prism changes the handedness of the incident beam carrying 
OAM and causes a sign change in its winding number. 

two-component superposition (a \£) +j3 |— £))/\/\a\ 2 + |j8| 2 . 




FIG. 2: Scheme for creation of superposition of the OAM states: 

(\+e) + \-£))/V2. 

With the mathematical details discussed at length in Ap- 
pendix|A] the process can be described as follows. The initial 
state of the light field before entering the first beam splitter is 
\£). After the passage through the first beam splitter, the beam 
is equally split into two paths and the state is a \£) j + j3 \£) 2 - 
Beam 1 is the part of the beam that passes through the beam 
splitter and beam 2 is the one that is reflected. Beam 2 later 
passes through the dove prism to undergo the transforma- 
tion \£) — > |— I) as discussed above. Thus, just before en- 
tering the second beam splitter the state of the light field is 
a\£) l +j8 \-£) 2 , which becomes a \£) + j8 \-£) after the re- 
combination in the second beam splitter. The complete Mach- 
Zehnder interferometer can be described via the transforma- 
tion: 

fu \£)\ ( r 2 u \-£)-t 2 ^u \£) \ 

\ J \itruo\-e)+irte i t>u \£)J ' K ' 

With the choice of 50-50 beam splitters we have r — t — 1 / V2 
and <p = K to obtain the state at the output ports of the inter- 
ferometer 

1 u ( H) + \i) \ (CS 

Thus, by ignoring the port 2 and using the state from port 1 
we obtain the required superposition state. 

To generate arbitrary superposition of the kind (a \£) + 
/3 y Ot 2 + /3 2 scheme in Fig. [2] needs to be changed 
slightly: (i) the input state needs to be a normal Gaussian 
beam with no orbital angular momentum, (i) the dove prism 
needs to be replaced by a hologram to transfer the gaussian 
beam into the state and (iii) an extra hologram is added 
into the lower path to obtain conversion from the Gaussian 
beam into the state \£). 



A superposition of two OAM states can be detected via a 
photo-detection scheme by taking a look at the spatial profile 
of the characteristic interference of the beams that are part of 
the superposition. For example, the interference pattern cre- 
ated by the superposition of \l) and |— £) is comprised of 21 
bright lobes that are equally distributed along the azimuthal 
direction. The interference pattern of arbitrary superposition 
would contain \£\\ + \£ 2 \ lobes where l\ and £ 2 are the com- 
ponent OAM states. Thus the observation of the described in- 
terference pattern guarantees the creation of the coherent su- 
perposition of the appropriate states. Further details offered 
in Sec. [IV] for detection of BEC vortex superposition can be 
applied here as well as the phase structures of the OAM super- 
position in light and vortex superposition in BEC are identical. 



III. COUPLING OF OPTICAL VORTEX BEAMS TO BEC 

In this section we discuss two methods to transfer the gener- 
ated superposition of optical orbital angular momentum states 
to the states of atoms. Transfer of pure OAM states of light 
to BEC has been studied by Marzlin et al. [17] and Nandi 
et al. ifToh . Nevertheless, for the transfer of OAM superposi- 
tions one needs a special transfer scheme. The internal level 
scheme of the atoms and the transitions of interest are shown 
in the Fig. [3] The internal hyperfine quantum number of the 
initial non-rotating state |0) is mp = and that of the final 
states |+) and |— ) is mp = 2. Thus, the internal quantum 
number of both the final states is chosen to be the same to 
obtain a pure superposition of the vortex states while all rest 
of the quantum numbers needed to describe the two states are 
exactly identical. The intermediate states, |/) and \i'), through 
which Raman coupling between the initial and final states oc- 
curs, have the hyperfine quantum number mp = 1. The two 
components of the optical OAM states correspond to the Rabi 
frequencies £2 + and £2 and both are 0+ polarized. The cou- 
pling field is designated by the Rabi frequency £2 C and it is a 
polarized. The OAM of the state |0) is zero and that of |/) and 
|+) is +£ and that of |/') and |— ) is —£, where £ is a positive 
integer. Even though we are targeting generation of superpo- 
sition of two counterrotating components, it can be quickly 
seen that this scheme can also be used as is for transfer of ar- 
bitrary superposition of two arbitrary OAM states of light to 
the BEC cloud. 

Due to the various possibilities available for the internal 
state quantum numbers and the OAM quantum numbers there 
are five distinct states that the atoms in the BEC can have as 
a result of their interaction with the optical fields. Thus, the 
BEC cloud can be described through a five-component spinor 
{^o^+^-^i^ji}. All the internal states are assumed to 
be trapped by the trapping potential and thus there are mu- 
tual interactions among all the components. The optical fields 
couple various components of the spinor BEC, as shown in 
Fig. [3] Thus, we can write the evolution equations for the 
spinor components 

i * = £2;^ + £T ^ + (JG/K) ^0 
i*+ =£ll^, + {,m/Ti)^ + 
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FIG. 3: Internal level scheme of the trapped BEC atoms and the ap- 
plied fields. The fields are specified through the corresponding Rabi 
frequencies and the polarization state. The fields are highly detuned 
from the excited states to avoid populating the excited states and en- 
sure the stability of the condensate. The two-photon detunings S 
are needed for one method of transfer of population from the non- 
rotating state to the vortex states but needs to be zero for the STIRAP 
like transfer mechanism. 

i* = Q*W f + (J^/h) } if- 
i = -AWi + £l c W+ + il+^o + (JG/ti) % 
i % = — A 1 ?,/ + Q c ¥_ + Q_¥ + (Jfi/h) *S t ■ ( 6 ) 

It is important to note the time derivatives appearing on the 
left-hand side of the equations. In general, the spinor com- 
ponents have both spatial and time dependence; nevertheless, 
for the optical evolution the time dependence is of prime im- 
portance. The optical coupling is governed by the Rabi fre- 
quencies Q. c (for the coupling field) and D.± (for the OAM su- 
perposition). Here, A is the single-photon detuning. In setting 
up the above equations we have assumed that the two fields 
are in two photon resonance with the transitions |0) — |+) and 
|0) — |— ), meaning the parameter 8 appearing in Fig.|3]is zero. 
This detuning can be added later on, if needed, as will become 
clear later. Furthermore, 

^ = (^ + r- A i) + 7? (|^ | 2 + |^ + | 2 + |vp_| 2 + |vp ! f + |^,,| 2 ) 

(7) 



where is the Kinetic energy operator, "V is the trapping 
potential of the BEC and \i is the chemical potential. The 
interatomic interaction strength is given by r\ = 4nh 2 aN/m = 
NU, where a is the s-wave scattering length, N is the mean 
number of atoms in the cloud, and m is the atomic mass of the 
atom. The absolute square terms in Eq.(|7]l thus correspond to 
the self-energy arising due to the interparticle interaction. 

Noting the spatial profile of the LG beams, we obtain 



a ± (r)=a ± a c- r2 / w2 ^ e±We ife (8) 



with £Iq being the Rabi frequency of the atom-field interac- 
tion. a± are the probability amplitudes of the vortices in 
the superposition state we want to create. We assume that the 
Gaussian fall-off of the light intensity is over a length larger 
than the transverse size of the BEC cloud. As a result, we can 
ignore the gaussian term in the LG beam profile to arrive at 

n ± = fl± (V2)i f in (^yV. (9) 

Note that we have incorporated the phase term inside the com- 
plex quantity x±iy = rcos0 ±isin0 = re ±u ^. Similarly the 
gaussian fall-off of the intensity of the coupling beam (or D. c ) 
along the transverse direction is ignored as the BEC cloud size 
is much smaller than the transverse beam profile. 

Now, we perform complete adiabatic elimination of the ex- 
cited levels \i) and |z'). We are completely justified in using 
the adiabatic elimination procedure as: (i) the optical evolu- 
tion happens on a time scale much faster than other processes 
in the BEC; (ii) the optical fields are so detuned from the ex- 
cited states that the excited states are never occupied. Thus, 
by substituting ^P, = 4V = and eliminating ^ and Sfy from 
the Eq. ©we arrive at 



i^ 



i*o = f(^ + r- A i)^o + ?(| l Po| 2 - 



1^4 



|¥_| 2 )¥ + -(|n 4 



A 



i(^ + r-M)^+ + y(^+l 2 + l*-l 2 + l^ol 2 )^+ + 7iaf l P+ + 7n + £i f * i Po 

n n A A 



i* = i + y - n)y_ + J (|*p+ 1 2 - 

n n 



(10a) 
(10b) 
(10c) 



Our ultimate interest lies in the temporal dynamics of the 
populations of various components of the BEC. Thus, the spa- 
tial part of the condensate is to be integrated out of the above 
equations. To deal with the spatial profiles of the spinor com- 
ponents of the BEC explicitly, we note that the general state 
of the BEC cloud can be written as 



where the spatio-temporal projections are given by 

*F (r,f) = a(t) (r|0) = a(r)exp[i Qi/h-K)t] yr g (r) 

(r, t) = p (t) (r|+> = p (t) exp[i (5 + p/h - k)t] r) 

¥_(r,t) = 7(0 (r|-> = tft)e*p[i(8 + n/h-K)t] \^-(-1,t). 

(12) 



(r|¥) - a(t) (r\0)+P(t) (r|+) + y(r) (r|-) (11) The two . p hoton detuning 8 introduced in the phase factors is 
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defined as 8 — V± — V c — (Oo±, where v± and V c are the an- 
gular frequencies of the optical field carrying OAM and the 
coupling field respectively, and ©o± is the energy level differ- 
ence between the states |0) and |±) expressed as an angular 
frequency. Furthermore, the parameter K, in general, depends 
on the interaction strength Tj between the BEC atoms and the 
dimensional parameters of the BEC cloud, which in tern de- 
pend on the trapping potential. For example, for the Harmonic 
trapping potential K is given by 

nha sc N _ 77 
m{2nyi 2 h L\L Z 4(2n)y 2 h L\L Z ' V ; 

It can, however, be noted that as we are dealing with conden- 



sates consisting of single species, K need not be explicitely 
incorporated in the phase factors as it just causes a shift in the 
energy equally for all the component states of the BEC cloud 
and does not affect the population dynamics in a non-trivial 
manner. 

The ansatz for the spatial profiles of the non-rotating ground 
state and the vortex states depend on the type of trapping po- 
tential we choose. Before discussing different trapping po- 
tentials we take our general formalism a little further. After 
substituting the ansatz (fT~2T > into the equations ( fTOb we obtain 
the temporal dynamics of the populations of the BEC compo- 
nents of interest: 



i cc(t) - (| - k) a(tj\ y/ g (r) = -(3- + r-n) a (t) % (r) 



^(|o(0| 2 |%(r)| z +|j8(/i|-|^ ('■.r.l| J + |y(/)| J </A l-M-.hal/.y/.lr! 



\flr 



I 2 ' 

11 



w 



+ ^«e- 1 ^j3(f)rv + (^r)+fl*_e^ 7 (f)¥v-M,r)] I \ %£l c , 



w 



ij8(0-(|-»c+5)j8(0]^v + (^r) = i(^ + r-M)j8(r)^v + (^r) 



|a( f )|Vg(r)| 2 + |j3(0lVv + (^r)| 2 + |7(0lVv-(-^r)| 2 iM*,r)/3(f) 



■~\a c \ 2 P(t)y y+ (W + ^ + e U,> (^) Oo^aCOVkW. 



i7(0-(f -Jc+5)7(0]^v-(-^r) = i(^ + r-M)7(0^v-M,r) 



|a(0rl%Wr+l/5(0rlv^+(^r)r+|y(0lV-(-^')rJ^-(-^i-)y(0 

I 

(14) 



i|£2 c | 2 7(0rv-M,r) + Ve-^ £l £l* c a(t) % (r) . 



Note that the exponential factors exp[i(/x//i — K)t\ and 
exp[i(5 + n/Ti — K)t\ have canceled as expected. Also note 
that I is taken to be a positive integer and ± signs are used 
explicitly to easily distinguish the two counter-rotating vor- 
tex components; we will follow this convention from here on- 



wards to avoid excessive use of the modulus sign. Now we 
formally perform the coordinate integrals to remove the spa- 
tial dependence. In terms of the resulting integrals, which are 
summarized in Appendix B, the rate equations are 



ia(f) = (r, + V>(f) + (/,>(^ 
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i$(t) = [T + +V + (t) + 8}IS(t) + (I g+ (£) \a(t)\ 2 +I ++ (£) \P(t)\ 2 + I + 4£) |y(r)| 2 )/3(f) 



-m- 



iy(t) = [T_ + V_{i) + 8]y(t) + (I^{l)\a{t)\ 2 + I + ^)\m\ 2 + (15) 



r 



Note that the spatial integrals are denoted by the letter / with 
different subscripts, superscripts and in most cases with the 
argument of I to properly distinguish them. These integrals 
are explicitly written out in the Appendix B. 

To note, the above set of equations is very general and ap- 
plicable to wide variety of trapping potentials, provided the in- 
tegrals are evaluated by appropriately taking the spatial forms 
of the wavefunctions, W g (r), y/,,±(r), suitable for the trapping 
potential. Furthermore the equations can also be applied to 
situations where the applied OAM and coupling fields have 
a certain time profile. In case the spatial ansatz for the BEC 
wavefunction in the Thomas-Fermi limit then the kinetic en- 
ergy terms can be ignored as we do in the case of the mexican 
hat potential. The application of these equations shall become 
clear as we discuss various cases below. 

In the forthcoming subsection we look at two different 
trapping potentials and two different transfer mechanisms of 
OAM to a BEC. 



A. Harmonic Potential Trap 

The harmonic potential we consider is of the form, 



1 



We assume a pancake shaped BEC cloud for which co± < C0 Z 
and as a result, L± > L z . The spatial wavefunctions for the 
BEC trapped in this kind of trap can be taken to be, 



V/,(r)=exp{-(l/2)[(r/L ± ) 2 + (z/L z ) 2 ]}/(^ /4 ^ 1/2 ) 



: r ri e 



%to/(vWf') 



(17) 



V(x,y,z) = -m[(o{(x l +/) + CO z z 



(16) 



The treatment, so far, is very general and would work for any 
value £ of the OAM of the incident light. To understand the 
dynamics more clearly we restrict ourselves to a particular 
value of the OAM, £ = 2. However, one may note that even 
though the spatial integrals would have different values for 
different £ values, the general idea remains the same. 

We start with Eq. $15[ and substitute the spatial integrals for 
the particular case of the 3D harmonic trap, which are listed 
in Appendix B, to arrive at: 



ia(t) = -K-a(f) + (ioj z + io) ± )a(f) + ia) ± a(f) + 3K-|a(f)| 2 a(f) + K-(|a(f)| 2 + |j3(f)| 2 + |7(f)| 2 )a(f) 



T °o. , 
A \m(0±w 



- ' A m(i>\w L 



i/3(0 H-S - k))3 (r) + + ^ ± )/3 M ^ ^ 
A A mOJiw 

ir(0 - (-5- ^)r(0 + (^«z + |co^)7(0 + |^7(0 + ^(l«(0l 2 + |y3(0l 2 + l7(0l 2 )7(0 + |^(l)8(0l 2 + l7(0l 2 )7(0 



i|afy(0 + ^oa:«- 



mCO±w 



■a(t). 



(18) 



Note the appearance of the Harmonic oscillator trap parame- 
ters. At this stage, we note that, |a(f)| 2 + |/3(f)| 2 + |y(0l 2 = 1> 
and eliminate terms that are common in all the equations as 
they do not give any non-trivial contributions to the popula- 
tion dynamics of the system. We also add an extra assump- 



tion to relate the properties of the BEC trap and the profile of 
the OAM carrying light beam such that 2y/2h/(m(Oj_w 2 ) = 1. 
Note that this only scales the optical coupling constant ap- 
pearing in the equations to a manageable number and does 
not change the physics in any way. Thus, we arrive at a gen- 
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eral set of equations coupling a non-rotating state of BEC to 
two counter-rotating vortex states carrying charges 2 and —2 



ia(r) 
i/3(0 
iy(t) 

Starting with these equations, we develop two methods for 
transfer of optical OAM to a BEC. It needs to be noted that 
the dynamics, as seen in the above equations, cannot have a 
steady state due to large detunings from the excited states and 
exclusion of radiative decays. This regime of parameters is, 
of course, necessary to make sure that the BEC atoms do not 
populate the excited state, where the BEC may not survive. 
Thus, to obtain the complete population transfer in a robust 
manner we need to construct dynamics that would not cause 
the population to oscillate between the initial state and the fi- 
nal states but would rather show a one-way trend in the direc- 
tion the population is moving. This, of course, could be done 
in both a time-reversal symmetric or non-symmetric manner. 
We study both these avenues. The non time-reversal symmet- 
ric scheme involves linear chirp of the coupling field, meaning 
a time dependent two-photon detuning 8. The time-reversal 
symmetric scheme is based on the traditional STIRAP popu- 
lation transfer mechanism. To note, the final state is not a sin- 
gle state but a superpositon of two states. We consider these 
two techniques one by one in the following subsections. 



1. Superposition of Vortices in BEC via frequency chirp 

To study the generation of the superposition of two counter- 
rotating vortices in BEC we start with Eq. (fT~9b and introduce 
two simplifications given by Q.q = Q. c and |£2o| 2 /A = a>±. 
There is no loss of generality in these assumptions and appear- 
ance of any other constant that related the light beam intensity 
parameter with the transverse trap frequency would demon- 
strate the same physics except the time scales may be a little 
different. Our aim here is, of course, to demonstrate that su- 
perposition of vortex states can be created in a robust man- 
ner despite the presence of inter-particle interactions as nor- 
mally observed in BEC. The same formalism can be applied 
to the experimental situations to rigorously determine the time 
scales over which the population transfer occurs. In any case, 
the optical time evolution time scales (/xs) are much shorter 
than the time scales for spatial evolution of the BEC (ms). 

We obtain a simple set of equations governing the popula- 
tion of the three components of the BEC: 

ia(t) = 3K\a(t)\ 2 a(t) + (o_ L (a* + P(t) + a*iY(t)) 



respectively: 



(19) 



i$(t) = (S + 2co ± )l3(t) + \K(\P(t)\ 2 + \7(t)\ 2 )m 
+a)±a + a(t), 

iy(t) = (5 + 2a) ± )7(0 + ^(l^(0l 2 + l7(0l 2 )7(0 

+(Oj_a-a(t). (20) 

The optical vortex superposition can be transferred to the BEC 
via a continuously chirped control pulse Q. c . The chirp is mod- 
eled by introducing linear time dependence in the two-photon 
detuning in the form 8(t) = C(l — Clot), where C is some ap- 
propriate constant. We discuss the solution of the above equa- 
tion in Fig. [4] 



1.0 



k*)P \\ 

; Mf)l 2 U, 














0.0 0.5 1.0 15 2.0 



fi * 

FIG. 4: Generation of the vortex state superposition via the linear 
frequency chirp technique. An equal superposition of the component 
vortex states is generated here starting with a non-rotating state of 
the BEC. The time plotted on the x-axis is measured in seconds. The 
linear chirp used for complete population transfer is given by the 
two-photon detuning 8(t) = C(l — Slot), where C = 2£Iq and flrj = 
3kHz. Other parameters are a>x = 132 Hz, a = 5 nm, L ± = 2.35/xm, 
L- = lAfim and K = 1.7kHz. 

The solid line in the plot is the transfer function 

F(t) = \a{t)\ 2 -\fl{t)\ 2 -W)\ 2 (2D 

which varies from the value 1, when all atoms are in state |0), 
to the value —1, when all the atoms are in state ( \/3 | +2) + 
\/2\— 2))/v / 5, which is the required vortex superposition for 



= 3K\a(t )\ 2 a(t) + i|n | 2 a(f) + ~Q5Q c «j3 (t) + a*_ y(t)) , 
= (8 + 2(o^(t) + l -Km)\ 2 + \7{t)\ 2 m) + \\&c\ 2 Ht) + \p. a* c a + a{t) , 
= (S+2a> ± ) Y (t) + ~K(|j3 (r )| 2 + \y{t)\ 2 )y{t) + \\a c \ 2 y(t) + ^ Q £l* c a_a(t) . 
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60:40 division of the population among the two vortex com- 
ponents. The population of state |0) is shown by a dashed 
line and it varies from 1 to 0; whereas the populations of the 
states |+2) and |— 2) denoted by |/3(f)| 2 and |y(f)| 2 are given 
by dotted line which vary from to 0.6 and 0.4 respectively, as 
expected. The coupling field frequency is to be varied linearly 
such that it sweeps from one side of the two-photon resonance 
to the other side to facilitate complete population transfer. For 
modeling purposes we vary the two-photon detuning 8 so that 
the system sweeps through the two-photon resonance. 

An alternative method we discuss below, to accomplish the 
OAM transfer, is similar to the traditional STIRAP technique 
of the counter-intuitive pulse sequence. This technique has 
been applied for the transfer of a pure OAM state to the BEC. 

2. BEC vortex superposition via STIRAP-like pulse sequence 

The only change needed here from the previous subsection 
is to make the Rabi frequencies time dependent and to deploy 



counter-intuitive pulse sequence that is required to obtain the 
population transfer. 

We add the time dependence to the Rabi frequencies by via 
Qo = |flo|/(0 and a c =\Oo\g(t), (22) 

where the temporal profiles of the beams are taken to be of the 
form: 

/(f)=/oe~(^r) , g(r)=g e~(^) . (23) 

With these modifications and that 8 — the rate equations (T% 
become 



ia(0 = ^|n | 2 /(0 2 «(f) + 3 J c|a(0| 2 a(f) + l|£lo| 2 /(f)^) «/3(r) + a*_y{t)) , 

ij8( f ) = 2^y3(0 + I,c(|/3(0| 2 + l7(0l 2 ) y3(0 + l|^ | 2 g(0 2 y3(0 + 1|^ |V(0^(0«+ «(r) , 

iy(0 = 2«x 7(0 + ^|^o| 2 ^(0 2 7(0 + ^rc (|y3(r)| 2 + |r(0l 2 ) 7(0 + ^|^o| 2 /(05(0«- . (24) 



These equations can be solved in a straightforward man- 
ner. The typical results are summarized in Fig. [5] It can be 
seen that the transfer function plot signifies a complete trans- 
fer from initial state |0) to a 60:40 superposition of the vor- 
tex states |+2) and |— 2). Any combination of superpositions 
(50/50, 60/40 or 80/20) can be generated via the counter in- 
tuitive pulse sequence — the coupling-field pulse (time profile 
g(t) and pulse center ti) comes before the OAM superposition 
pulse (time profile /(f) and pulse center fi). 

To understand the pulse overlap necessary to obtain com- 
plete population transfer from the initial non-rotating state to 
the final vortex superposition state via the STIRAP process, 
we plot the transfer function F(t) vs the distance between the 
pulse centers t\ — 1% in Fig. [6] The value of F(t) close to 1 
means the population is primarily in the non-rotating ground 
state and the transfer is not efficient. Whereas, F(t) = — 1 
means complete population transfer to the vortex superposi- 
tion has occurred. The pulse parameters are chosen to be ex- 
actly identical; therefore, the distance between the pulse cen- 
ter could be taken as the measure of the pulse overlap. For 
the light pulses arriving at about the same time the transfer 
is inefficient. However, the population transfer improves as 
the pulse separation increases a little. One obtains complete 
transfer for the range {0.3,0.5} for the separation between the 
pulses; beyond that the pulses are so far away from each other 



that the STIRAP process does not work. 

We would further like to point out that our numerical study 
suggests that within the parameter range used to obtain Figs. [4] 
and [5] the excited states populations are negligible and we 
are justified in using the adiabatic elimination of the excited 
states. 



B. Mexican hat trapping potential 

We now study the transfer of OAM to a BEC cloud when 
the trapping potential is shaped like a Mexican hat (Sombrero) 
in the x — y plane (See Fig. [7]). The advantage of this poten- 
tial is the toroidal symmetry it offers. The toroidal trap con- 
figurations naturally support toroid shaped BEC cloud and as 
a result can sustain large vortices without disintegration into 
several single charge vortices. This offers stability to the vor- 
tex superpositions we are aiming to generate. The trapping 
potential, in this case, is of the form 13111 

V(p,z) = --Gmcolp 2 + -X f — ^ J p 4 + -mco 2 z 2 (25) 

where a and X are dimensionless parameters. The potential 
is harmonic in the z-direction. The dynamics of the OAM 
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-0.5- 



FIG. 5: Generation of the vortex state superposition: Results of the 
numerical solutions of the equations for STIRAP scheme shows the 
superposition 60:40 of the |+) and |— ) vortex states. The pulse pro- 
files f(t) and g(t) are of the form in Eq. i23\ with <y\ = 02 = 0.25 
and t\ = 1.0 and ti = 0.5 in the units of 1/£2q with fo/go = 0.5 to 
obtain the complete transfer of the non-rotating BEC into the vortex 
superposition. (See Eq. l!24t.) Other parameters are Qq = 2 x 10 5 Hz, 
A = lOfio an d 5 = 0. The quantity on the x axis is the scaled time 

n f. 



c 
o 

'*-> 

u 
c 
3 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
Distance between pulse centers (scaled time) 

FIG. 6: Plot of the transfer function (Eq. J2U ) versus the distance 
between the pulse centers t\ — £2. Note that t\ — ti > (or t\ > £2) 
suggests the counter-intuitive pulse sequence needed for the STIRAP 
based population transfer process. The other parameters are identical 
to those noted in the caption of the Fig. [5] 



transfer to BEC can be studied in the Thomas-Fermi (TF) ap- 
proximation where the kinetic energy of the BEC cloud is ne- 
glected. In this approximation, the spatial part of the BEC 
wavefunction is given by 



Max 
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/ 00005 



FIG. 7: The mexican hat potential in two dimensions with a = 2.0 
and X = 0.005. 



Here we have assumed that the wavefunction in z-direction 
has a Gaussian form and has little effect on the vortex dy- 
namics of interest in the transverse direction [18]. The spatial 
profile of the above wavefunction is shaped like a toroid (or a 
donut) with a hole in the center. Thus, the two radii (inner and 
outer) are necessary to describe the shape of the BEC cloud. 
The radii can be found by setting | \jf(l,p, 0,z)| 2 = to deter- 
mine where the particle density goes to zero. The two real 
solutions for the radii are given by 



Rz 



oM±yJo 2 M 2 +4rlgii 



(27) 



where M — m(o\ and g = m (O^/h . A typical spatial form 
of the ground state BEC cloud in the mexican hat potential is 
shown in Fig. [8] along one transverse direction (i or y). The 
toroid shaped particle distribution should be clear from the 
figure. 




FIG. 8: The Thomas-Fermi spatial wavefunction of the BEC in a 
Mexican Hat potential in one dimension. 

The spatial part of the BEC wavefunction for different com- 
ponents is taken to be 



\j/ g (r) = y/(0,p,</>,z), 
Vv±(r) = \ir(±£,p,Q,z). 



(28) 



Using the above one can evaluate the spatial integrals appear- 
ing in the Eq. ( fT3T > governing the population dynamics. We, 
furthermore, ignore kinetic energy of the system as we are 
working in the TF approximation, in this case and choose the 
optical field frequencies so that there is two-photon resonance, 
i.e. the two-photon detuning 5=0. The light fields Q. c and 
Cl± have the same saptio-temporal form as discussed in sec. 
1111 A 21 Numerically solving the evolution equations, we get 
the transfer curves for the mexican hat potential depicted in 
Fig. [9] We observe that the time-scale of transfer is about 
the same as in the case of harmonic trap. The figure depicts 
generation of a 60:40 superposition of |+2) and |— 2) vortex 
states. Furthermore, it is important to note that any arbitrary 
superposition could be obtained in the mexican-hat trap; the 
generated vortex state would be naturally stable and would not 
disintegrate into single charge vortices. 




0.0 0.5 1.0 15 2.0 



n t 

FIG. 9: Generation of the vortex state superposition in a Mexican 
hat trap: Results of the numerical solutions of the equations for STI- 
RAP scheme shows the superposition 60:40 of the | +) and | — ) vortex 
states. The Rabi frequency £lo = 1kHz and A = 100f2(). 



IV. DETECTION OF THE BEC VORTEX SUPERPOSITION 

The resulting vortex-superposition state could be detected 
by imaging its particle density distribution, which is propor- 
tional to an interference of its components. Fig. \W\ shows 
the x-y cross section of this interference pattern for a partic- 
ular vortex state |¥) v = a \£ = +3) + /3 e ie \-l = -3). The 
spatial profile of the interference pattern can be obtained by 
evaluating | (r|W v ) | 2 = A[l + 2aj5 cos(2^ - 9)} where A = 
| x P v -i-(r, (j>) | 2 is just the toroid-shaped particle density distribu- 
tion of the vortex states. The "cosine" term dictates the mod- 
ulation of the particle density distribution along the azimuth 
with 2£ oscillations as the azimuthal coordinate <p changes 
from to 2n. The relative phase of the two components, 9, 
just causes offset( or rotation) of the interference pattern by an 
angle 9/(2£). Note that r is a two-dimensional vector signi- 
fying the position of a point in polar coordinate system {r, (j)}. 
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FIG. 10: Detection of the vortex-superposition through characteristic 
interference of a general normalized state a \ +£) + 13 e 1 e | — £}. £ = 3 
and 0=0 unless specified otherwise, (a) a 2 : j3 2 = 1 : 1 ;V — 1, 
giving m = 2£ lobes in the interference pattern, (b) a 2 : /3 2 = 0.1 : 

0. 9 ; V = 0.6. (c) a 2 : /3 2 = 0.9 : 0. 1 ; V = 0.6. Patterns (b) and (c) 
are identical, (d) Phase determination: a 2 : jS 2 = 1 : 1 , 8 = n\ V = 1, 
notice rotation of the pattern with respect to that of (a) by an angle 
0/m. After adding one more unit of OAM in the superposition to 
arrive at a \+£+ 1) + [3 e l6 |— £ + 1) distinct patterns (e) and (/) are 
obtained instead of (b) and (c). 

The visibility of such a pattern can be readily arrived at via: 
y = W-4™ =2aj3 (29) 

Imax hnin 

as lmax = A ( 1 + 2afi ) and /,„,„ = A ( 1 — 2afi ) are the extremal 
intensities. 

In general, for such a superposition the spatial profile of 
the interference pattern contains m = 2£ lobes. The visibility 
V = 2aj3, as determined above, gives a measure of the asym- 
metry in the amplitudes. The pair {(X,f}} can be determined 
using measured V and the normalization condition a 2 + f} 2 = 

1 . However, this does not assign the amplitudes to the states 
|+) or |— ) with certainty, as the patterns in (b) and (c) are 
identical. We propose shining a OAM £ = +1 light of cr + po- 
larization to obtain the vortex state a |3 + 1) + j3e l9 |— 3 + 1), 
the resulting interference pattern is shown in Fig. [10] (e) and 
(/), which now clearly differentiates between the two ampli- 
tude values that gave same visibility in (b) and (c). The phase 
difference 9 causes rotation of the whole pattern by an amount 
9/m as shown in (d) of the figure. Existing schemes for de- 
tecting vortex states j32tl could also be extended to detect a 
superposition of vortex states. 

V. CONCLUSION 

To summarize, we have studied two different mechanisms 
for transfer of superposition of optical angular momentum of 
light to vortices in BEC. An interferometric scheme for gener- 
ation of an arbitrary superposition of two different OAM states 
of light is also presented in great detail. We also discussed a 
couple of trap configurations and showed that the OAM trans- 
fer from light to BEC works independently of the trapping 
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potential and despite inter-atomic interactions present in the 
BEC. The applications of the techniques discussed here are to 
a memory for the OAM states of light with further applica- 
tions in the quantum or even classical communications proto- 
cols using OAM states of light. The superposition states of 
vortices in BEC could also be used as a qubit for quantum in- 
formation processing; however, further work is clearly neces- 
sary on that front. Further, various advantages of using atoms 
versus photons for interferometric metrology are discussed in 
Ref. 113 311 . The superposition of counter-rotating currents in 
the BEC could also be used for inertial sensing, especially 
as a gyroscope. Several advantages offered by a gyroscope 
based on superpositions of counter-rotating vortex structures 
would be the tunability of the effective de-Broglie wavelength 
by choosing appropriate atomic masses and the angular ve- 
locity and phase sensitivity via the choice of the quantized 
angular momentum of the atoms. 



written as 

Were r and t are the reflection and transmission amplitudes for 
the input at port 1, and r 1 and t' are the parameters for the port 
2 (See Fig. [TTJ. Noting that the matrix R should be unitary, 
meaning, we obtain 

{-,-:) <«> 

The determinant of a unitary matrix has a modulous of one, 
therefore we have rr 1 — tt' = However, this factor multi- 
plies all the elements of the matrix, thus we can safely choose 
it to be 1, given by the choice of y = 0. By equating the corre- 
sponding elements on the RHS and LHS of the above equation 
we obtain 
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APPENDIX A: DETAILED DESCRIPTION OF THE MACH 
ZENDER INTERFEROMETER 

In this appendix we offer mathematical details of the inter- 
ferometric scheme for generation of the OAM superposition 
for light. We suggest references lHH HH Hfl to an interested 
reader for further details on treatment of the optical elements. 




FIG. 1 1 : A schematic representation of an ideal beam splitter with 
input ports 1 and 2 and output ports 3 and 4. The amplitudes for 
reflection and transmission from port 1 are, respectively r and t. The 
corresponding amplitudes for a photon entering port 2 are r' and t' . 

We chose a basis state representation as 




(Al) 



In this representation the general beam splitter matrix can be 



r = r* and t' = -t* (A4) 

If we rewrite these factors as complex exponentials, |r|e l5r , 
jfje 1 *, |r'|e l5 '-' and |f'|e l5 ''. by division of the above equalities 
we obtain 




(A5) 



From dA4l ) we can see that \t\ = \t'\ and \r\ = |r'|; using this in 
the above equation we obtain, 

8, - 8 r + 8,i -8/ = n. (A6) 

For the case of a symmetric beam splitter, which has the same 
effect on a beam incident through port labeled 1 as on a beam 
incident through port labeled 2, r — r* and t = t ', and we have 
8, — 8, + 8,i — 8 r i = tc/2. Thus the transmitted wave leads 
the reflected wave in phase by n/2 radians. This is a general 
property of symmetric beam splitter ll34ll . This also implies 
that r and r 1 are purely real and equal and t and ?' are purely 
imaginary and equal to each other. Let r = r 1 = r and t = t' = 
i t be the Thus the beam splitter matrix becomes 

*-(£»). (A7, 

For the case of a 50-50 beam splitter, for which r — r* and 
t = t' and also |r| = |r'| = \t\ = \t'\, we have 

£(50-50) = _L (1 ^ (A8) 

Now we consider the complete transformations for the Inter- 
ferometer considered in Fig. 12 
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(A9) 



Observing the output state, it can be easily seen that a general superposition could not be generated by this method. Instead we 
need only one beam-splitter to be imbalanced, meaning \r\ ^ \t\, and the other one to be 50:50. However, we can still choose 
both the beam splitters to be symmetric, i.e., r = r J and t = t'. Thus using Eq. (IA7t and Eq. dA8b we can attain the output state. 
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(A10) 



Thus, by choosing t = a + and f = a_, with the condition that f 2 + I 2 = 1, and e 1 
a- \—£))/V2, which is a general superposition state. 



= — 1 we obtain the state uo(a + \£) 



r 



APPENDIX B: THE SPATIAL INTEGRALS 



trapping potential, are given below. 



The various coordinate integrals needed to eliminate the 
spatial parts of the BEC spinors to arrive at the population 
evolution equations are given below. The definitions given 
below are independent of the trapping potential or the ansatz 
used for the spatial profile of the wavefunction. 



T±(£ 
V ± (£ 

I g ±(£ 



lHl(£ 



^(r)^ g (r) d 3 r, 
<(r)r%(r) d 3 r, 



1 

h 
l 

h 

l -J y* ± (±£,r)J7y v± (±£,r) d 3 r, 
jJ V ;±(±e,r)-r Wv± (±£,r) d 3 r, 
Jj|y/ g (r)| 2 |y/ g (r)| 2 d 3 r, 



|y/ v± (^r)| 2 |y/ g (r)| 2 d 3 r, 
|y/ v+ (^r)| 2 |y/ v+ (£,r)| 2 d 3 r, 



±J |v/v-(^r)| 2 |v/v-(^r)| 2 d 3 r, 
^W(^) Vg(r)d 3 r, 



'y/2r 



. j tf ± (£,r)^ y/ g (r)d 3 r. (Bl) 



The values of the above spatial integrals, for specific vortex 
states corresponding to the charge of £ = ±2 and the harmonic 



T 

JH) 



h ( 1 2 \ l l a 

^U + ^H^+2^ 



'gg — 



h f 1 

n 



i 



^ 2 ) = 4^U + ^J = 4^ + 2^=^ 2 )' 



= 4jc, 



4? (2) 



/g(2) 



J gg 



(27T)2»L Z L ± 2 

4(27T)2/i L Z L ± 2 

8(27t)h L Z L ± 2 2 
■ 2L 4 x =2(h/ma>x) 



^ = 4 J (2), 



= ^=7^(2) =7^(2), 



& H \£) = V2Ll = m H \£) 



(B2) 



Note also that at this stage we have used various properties of 
the quantum harmonic oscillator as summarized below: 



2mL 2 2 



mffli , i.e., L± 



h h 

, and L z = a , 

mC0± " V mCO- 



The corresponding integrals for a the mexican hat trap are 
evaluated numerically for the parameter values a = 2.0 and 
A = 0.005. Also, as we deal with the Thomas-Fermi ground 
states for the mexican hat trap the kinetic energy is ignored. 

This completes the description of the spatial integrals, 
which allow us to obtain the time-dependent population equa- 
tions dT31 >. Those could be solved as discussed in the main 
text of the paper to study the transfer techniques for the vortex 
superpositions from the light to the atom. 
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